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The influence of hydrogen adsorption 
on magnetic properties of Ni/Cu(001) surface 
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Ni/Cu(001) is known as a unique system showing the spin-reorientation transition from an in- 
plane to out-of-plane magnetization direction when the Ni-overlayer thickness is increased. We 
investigate different relaxed multilayer structures with a hydrogen adlayer using the full-potential 
linearized augmented plane-wave method. The relaxed geometries, determined by total energy and 
atomic force calculations, show that H-monolayer strongly influences the interlayer distance between 
the Ni-surface and sub-surface layers yielding the outward relaxation of Ni-layer at H/Ni interface. 
Furthermore, large decrease of local magnetic moments at the top surface area is found for the surface 
covered by H. The magneto-crystalline anisotropy energies calculated for fully relaxed H/Ni-films. 
The spin-reorientation transition critical thickness of 4 ML is found in good quantitative agreement 
with the experiment. 

PACS numbers: 75.70.Ak, 75.30.Pd, 75.30.Gw 
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I. INTRODUCTION 

Ultrathin ferromagnetic films grown on nonmagnetic 
substrate show peculiar magnetic behavior 0]. One of 
the unique phenomena which is observed in the ultra- 
thin Ni films on Cu(OOl) substrate is the spin direction 
reorientation transition (SRT) from an in-plane to out- 
of-plane magnetization direction when the Ni-overlayer 
thickness is increased |^ |^ Q ■ It is very important for 
the spintronic magnetic device applications and its mi- 
croscopic understanding attracted recently both experi- 
mental and theoretical 0, interest. The gas adsorp- 
tion on magnetic film provides the way to monitor the 
surface magnetic properties resulting in strong influence 
on the SRT critical thickness dc- For Ni/Cu(001) films, 
the hydrogen adsorption is observed to reduce dc by w 
4 monolayers (ML) from its value of « 11 ML for the 
vacuum/Ni/Cu-films 

The key quantity which drives the SRT in ultrathin 
magnetic films is the magneto-crystalline anisotropy en- 
ergy (MAE). It determines the preferred film magneti- 
zation orientation by minimizing the free energy of the 
system. For the free-standing Ni/Cu(001) films, the dc is 
found to be determined 0, l3 by competition of the uni- 
axial MAE due to the tetragonal distortion of Ni-films 
(so called "volume" MAE) and the surface MAE. The 
adsorption of a gas produces non-trivial changes in ge- 
ometrical structure as well as magnetic properties of Ni 
films 0, influencing both volume and surface MAE, 
and results in change of the SRT critical thickness. 

In recent years it has become possible to make use of 
ah initio density functional theory to predict the MAE in 
ultra-thin fllms 0- From theoretical and computational 
point of view, the MAE calculations are extremely diffi- 
cult due to the very high energy resolution in the range of 
few /ieV which is required. For the II/Ni/Cu(001) films. 



the problem becomes even more complex due to the need 
of accurate account of structural relaxation. It requires 
us to use the most accurate total energy full-potential lin- 
earized augmented plane wave method (FP-LAPW) 
in order to account simultaneously on equilibrium geo- 
metrical and magnetic structures. 



II. METHOD AND RESULTS OF 
CALCULATIONS 

The experiments show that H2 adsorbs dissociatively 
on Ni(OOl) surface in fourfold hollow sites There- 
fore, as a structural model for ab-initio calculations we 
use free-standing Ni-films (H/Nid/H, d = 1-11) and Ni- 
films on Cu-substratc (H/Nid/Cuy/Nid/H, d = 1-6) with 
ordered p{l x 1)H adsorbate overlaycr as shown in Fig. 
1. 
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Figure 1: Surface layers of H/Ni,2/Cu(100) system. 

The in-plane experimental lattice constant of Cu 
acu = 3.615 A pH is used which remains unchanged in 
the calculations. All interlayer spacings dij (see Fig. 1) 
are relaxed to the equilibrium values. In order to deter- 
mine the equilibrium slab geometry we employ the FP- 
LAPW method in FLAIR implementation (unpublished 
improved and rewritten version of the original FLAPW 
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Figure 2: Relaxation of interlayer distances for 
H/Ni„/Cu(001) for n = 3-6 and for n = oo. For 
clarity, the data sets are shifted by 0.1 a.u. Lines serve as 
guide to the eye. 



codes 0)- The scalar-relativistic atomic force technique 
is employed for the total energy minimization. Here, 
the Perdew and Wang jl^ approximation for exchange- 
correlation potential is used with the plane wave energy 
cutoff Ecut of 13 Ry, and the 28 special k-points in the 
1 / 8th irreducible part of 2-dimcnsional Brillouin zone 
(2D BZ) arc used for the BZ integrations. The con- 
vergence better than 1 x 10"^ e/(a.u.)'^ is achieved for 
charge/spin densities, and better than 0.2 mRy/a.u. for 
the atomic force acting on individual atom. 

The relaxed interlayer distances for different free- 
standing films (H/Ni„/Cu7/Ni„/H, n = 3-6) are shown 
in Fig. 2 and compared with the relaxation of tetragonal 
p(l X 1)H Ni(lOO) surface. Clearly, we can distinguish the 
outward relaxation of top Ni layer as well as the strong in- 
fluence of Ni/Cu interface on inter-layer distances in the 
close vicinity of the interface. The relaxation of Ni/Cu 
subsurface is different only for 5 ML Ni. The deeper in- 
terlayer distances (thicker film) approach the bulk value. 



Table L Interlayer distances in a.u. 





HNiiiH 


HNieCuTNieH 


di2 


0.63 


0.60 


d2'i 


3.30 


3.28 


d-iA 


3.18 


3.16 


d45 


3.20 


3.19 



The equilibrium values of inter-layer distances for two 
model systems: eleven-layer tetragonal Ni-film covered 
with H and seven-layer Cu-film with overlayer of six Ni 
layers covered with H are shown in Table 1. Both in- 
terfaces (H/Ni, Ni/Cu) influence strongly the interlayer 
distances in the overlayer slab. Note, that the calcu- 
lated bond length du-Ni of 3.48 a.u. is slightly shorter 
than the bond length for the bulk nickel hydride (3.52 
a.u.). The outward relaxation of top Ni-layer is found, 
the H-ML influences strongly only the interlayer distance 
between the Ni-surface and Ni-subsurface layer. Below 
the interface, the interlayer distances oscillate around the 
bulk value for strained tetragonal Ni bulk {d± = 3.20 
a.u.). The electron screening in the metal is responsible 
for fast damping of these oscillation. 



Table IL Magnetic moments in 





HNiiiH 


HNieCurNieH 




0.238 


0.239 


Nia 


0.593 


0.588 


Nis 


0.640 


0.636 




0.633 


0.633 



The layer-resolved spin magnetic moments il/g for 
these systems are shown in Table 2. These values cor- 
respond to the magnetic moments in the Ni "muffin-tin" 
spheres [Rmt = 2.2 a.u.). It is seen that there is strong 
reduction of the spin magnetization for the top Ni-layer 
due to the interaction with the H-adlaycr. The strong 
hybridization of the H s state with the filled majority 
Ni d band changes the band structure and the surface 
density of states (SDOS) leading to the decrease of spin- 
majority and increase of the spin-minority SDOS. Away 
from the interface, the spin magnetic moments are slowly 
converging to their bulk values. The layer-resolved mag- 
netic moments for different H/Ni„/Cu7-films (n = 3-6) 
on Cu-substrate are shown in Fig. 3. We note that for 
very thin H/Ni„/Cu7-films with n=l;2 the Ni-local mo- 
ments disappear and the system becomes non-magnetic. 
With the increase of the Ni-film thickness, the local Ni 
Ms of 0.24 //B is formed at the Ni-interface and then in- 
creases away from the H/Ni interface. When approach- 
ing the Ni/Cu interface, the Ni-atom magnetic moments 
start to decrease again and become ~ 0.45 for Ni/Cu 
interface layer (depending of the Ni-film thickness) . 

The anisotropic energy density of a tetragonal ferro- 
magnetic film is written as [l^ : 

E/V = -Klml ~ K^mt - A^m^^^ 

-^(X>,^ + K^mt + Kimlml), 

where, terms are 2nd and 4th-order volume-type 
anisotropy constants, K"^ terms are surface/interface- 
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Figure 3: Layer-resolved magnetic moments of 
H/Ni„/Cu(001) for n = 3-6 and for n = oo. For 
clarity, the data sets are shifted by 0.1(n — l)/^^. Layers are 
labeled S, S-1, ... , starting with the surface layer S. Lines 
serve as guide to the eye. 

type 2nd and 4th-order anisotropy constants, rux^y^z are 
magnetization cosines vifitli respect to the crystal axes, 
and d is a thickness of magnetic film. We assume that 
the 4th-ordcr terms in Eq.(l) arc significantly smaller 
than 2nd-order uniaxial anisotropy constant A'J" , and ne- 
glect them. The MAE then can be characterized with 
the difference in the total energy when magnetization is 
oriented along [100] (in-plane, ||) and [DDI] (out-of-plane, 
_L) axes (MAE = £;[1GG] - E[001\Y 

We use the relativistic version [lj| of the FP-LAPW 
method to solve self-consistently the Kohn-Sham-Dirac 
equations with spin-orbit coupling included to obtain the 
ground state charge and spin densities for the magne- 
tization directed along [GGl]-axis. The fully optimized 
interlayer distances as discussed above are used for free- 
standing Ni-films and Ni-films on Cu-substrate with H- 
adlayer. The MAE is obtained by app lying the force 
theorem to the spin - axis rotation [13: from the self- 
consistent ground state charge and spin density obtained 
for the [001] spin axis, a calculation of the band structure 
for [IGG] spin axis orientation is performed, and difference 
of the single particle eigenvalue sums is then taken to be 
the MAE. For the MAE calculations the k-points mesh 
equivalent to 6400 k-points in the full 2D BZ is used 
guaranteeing the MAE convergence better than IG fieV. 

The MAE as a function of the Ni-film thickness for 
the H/Nid/H films with d = 3,5,7,9 ML is shown in 
Fig. 4. For the c? = 3 we found small and negative 
MAE (-G.02 meV) keeping the Ni-film magnetization in 
[IGG] plane. It becomes positive for the d = 5 (G.432 
meV) resulting in the out-of-plane magnetization switch- 
ing. The MAE is positive with further increase of d and 
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Figure 4: The MAE (circles) and MAE-f-SAE (diamonds) for 
the fully relaxed H/Nid/H films as a function of the film thick- 
ness d (in ML). 



shows pronounced oscillations. The linear interpolation 
MAE = d ■ Ky + 2 ■ Kj yields the estimates for the 
"volume" Kv= 0.1G6 meV/atom and "interface" Ki=- 
0.127 meV/atom MAE contributions. The calculated 
"volume" MAE agrees well with AV= 0.G835 meV/atom 
calculated for free-standing unrelaxed Ni-films without 
H-adlayer and extrapolated to T=G K experimental 
value of G.072 meV/atom of Ref. The H/Ni "inter- 
face" Kj is calculated to be substantially smaller than 
vacuum/Ni "surface" Ks — -0.447 meV/atom @ for free- 
standing Ni-films and extrapolated to T=0 K experimen- 
tal value of -0.7 mcV/atom Q. We note that for the very 
thin Ni-films considered here, validity of the linear fit for 
the MAE and its separation into "volume" and "surface" 
contributions is not well justified due to the strong depen- 
dence of the Ni-film magnetic properties on the thickness 
of the film, especially in the presence of the H-adlayer. 

In order to estimate the dc critical thickness for SRT, 
we take into account the shape anisotropy energy (SAE) 
due to the magnetic dipole interaction, which provides 
additional in-plane anisotropy. This anisotropy is esti- 
mated using the relation SAE = — 27rM^ to the spin 
magnetization density M (in CGS units). It yields for 
the SAE the values of -0.G16 meV, -0.G32 meV, -0.055 
meV and -0.073 meV corresponding respectively to the 
H/Nic;/H films with d = 3,5,7,9 layers (which are consis- 
tent with experimentally derived SAE/Ni atom of -0.GG75 
meV 0). This additional negative SAE slightly shifts the 
MAE in Fig. 4 downwards yielding the = 4. 

Indeed, the H/Ni^/H model is a way too simple to de- 
scribe quantitatively the SRT in in Ni/Cu layers with 
H-adlayer. As it was already mentioned in Ref. |^ 
the Ni/Cu interface can play an important role. In order 
to evaluate the infiuence of Ni/Cu interface on SRT we 
perform the calculations for H/Nirf/Cus/Nid/H {d = 3,4) 
films, where the 5 Cu layers play a role of the substrate. 
The use of 5 ML of Cu layers instead of 7 ML allows to 
reduce computational effort without producing any sig- 
nificant impact on the magnetic properties and the MAE 
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of H/Ni/Cu-films, since they are originated from Ni mag- 
netic film and not the Cu non-magnetic substrate, and 
the relaxed Cu-interlayer distance below the Ni-interfacc 
is found to be very close to its bulk value of 3.427 a.u. 
For the d = 3 we found negative MAE of -0.192 meV 
and SAE of -0.036 meV keeping the Ni-film magnetiza- 
tion in [100] plane. Already for c? = 4 the MAE becomes 
positive 0.543 mcV while the SAE is small and negative - 
0.052 meV. Again as in the case of free-standing H/Nid/H 
films we get the SRT critical thickness dc = 4. 

To summarize, we found that hydrogen adsorption for 
the Ni films on Cu substrate yields the reduction of the 
SRT critical thickness. The calculated dc = 4 for both 
H/Nid/H and H/Nid/Cus/Nid/H films agrees well with 
the experimental dc of 7 ML 0|. We show that the 
MAE has strong and oscillatory dependence on the Ni- 
film thickness which deviates substantially from the lin- 



ear fit. We attribute the decrease of the dc due to the H 
adsorption to originate from strong reduction of the mag- 
nitude for the Ni "surface" MAE contribution. In turn, 
this is caused by strong decrease of the exchange splitting 
at the H/Ni interface due to the strong hybridization of 
the H s state with the Ni-bands. 



Acknowledgment 

The financial support was provided by the Academy of Sci- 
ences of the Czech Republic (Grant No. A1010214), CMS Vi- 
enna (GZ 45.504), and by the RTN project "Computational 
Magnetoelectronics" of the European Commission (HPRN- 
CT-2000-00143). 



[1] S. Bliigel: Phys. Rev. Lett. 68 (1992) 851 and references 
therein. 

[2] B. Schuiz and K. Baberschke; Phys. Rev. B 50 (1994) 
13467. 

[3] R. VoUmer et al: Phys. Rev. B 60 (1999) 6277. 

[4] S. van Dijken, R. VoUmer, B. Poelsema, J. Kirschner: J. 

Magn. Magn. Mater. 210 (2000) 316. 
[5] C. Uiberacker et al: Phys. Rev. Lett. 82 (1999) 1289. 
[6] A.B. Shick, Y.N. Gornostyrev, A.J. Freeman: Phys. Rev. 

B 60 (1999) 3029. 
[7] E. Wimmer, H. Krakauer, A.J. Freeman: Adv. Electron. 

Electron Phys. 65 (1985) 357. 
[8] M. Weinert and J.W. Davenport: Phys. Rev. Lett. 54 

(1985) 1547. 

[9] H.J.F. Jansen: Phys. Rev. B 59 (1999) 4699 and refer- 



ences therein. 

[10] I. Stensgaard and F. Jakobsen: Phys. Rev. Lett. 54 
(1985) 711. 

[11] W.B. Pearson: A handbook of lattice spacings and struc- 
tures of metals and alloys. Pergamon Press, Oxford, 1958. 

[12] J.P. Perdew and Y. Wang: Phys. Rev. B 45 (1992) 13244. 

[13] L. D. Landau, E. M. Lifshitz, and L. P. Pitaevskii: Elec- 
trodynamics of Continuous Media. Oxford, 1995, p. 138. 

[14] A. B. Shick, D. L. Novikov and A. J. Freeman: Phys. 
Rev. B 56 (1997) R14259. 

[15] I. V. Solovyev, P. H. Dederichs and L Mertig: Phys. Rev. 
B 52 (1995) 13419 and references therein. 

[16] M. Farle at al: Phys. Rev. B 55 (1997) 3708. 



